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Based on a general form of the effective vertex functions for the decays of P-wave charmonia XcJ, 
compact formulas for the angular distributions of the J/ip in the XcJ J/ip*1 decay and the muon 
in the subsequently cascade decay XcJ — > J/tpJ M + M~7 are derived. While the former are the 
same as those obtained in a different approach by using the angular momentum relation, the latter 
ones are new. Similar to the J/ip case, the XcJ polarization observables can be expressed in terms 
of rational functions of the yields spin density matrix elements of XcJ- A reweighting method is 
proposed to determine the production spin density matrix elements of the Xc2- Some generalized 
rotation-invariant relations for arbitrary integer-spin particles are also derived. Special cases are 
applied to the polarizations of Xci and Xc2- In particular, we discuss the possible impact of Xcj and 
ip' feed-down on the prompt J/ip polarization. As a phenomenological illustration, we present the 
angular distributions in terms of the spin density matrix elements of Xci and Xc2 production in four 
different frames at the Large Hadron Collider. 
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I. INTRODUCTION 

The polarization of heavy quarkonium in hadroproduc- 
tion e.g. at the Tevatron and the LHC has been a long- 
standing issue in heavy quarkonium physics[l]. Non- 
relativistic QCD (NRQCD)[2], a rigorous effective field 
theory that characterizes the nonrelativistic feature of 
heavy quarkonium, gives the description that a QQ pair 
may be formed in a color-octet (CO) state during the 
hard reaction at short-distances before its hadronization 
into a color-singlet (CS) physical quarkonium by radiat- 
ing soft gluons. Hence, J /ip, a cc bound state with quan- 
tum number J PC = 1 , is produced dominantly by 
gluon fragmentation into cc[ S[ ] first and it then evolves 
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into S[ in the high transverse momentum region [3]. The 
gluon fragmentation mechanism guarantees that the J/ip 
is produced in the transverse polarization in its hclicity 
frame (HX) when its pt is sufficient large, considering 
that spin flipping is order of v 2 suppressed according to 
the power counting rule in NRQCD, where v is the rel- 
ative velocity between c and c. However, the data mea- 
sured by the CDF [4, 5] Collaboration at the Tevatron 
indicate that the J/ip is mainly unpolarized and even 
slightly longitudinal polarized at large px, say, up to 20 
GeV. The situation for ip(2S) is similar to J/ip. This is 
the polarization puzzle of heavy quarkonium production. 

The importance of charmonium polarization is also 
related to the efficiency uncertainty that the detec- 
tor acceptance for lepton pairs from the J/ip (or 
other heavy quarkonium) decay strongly depends on its 
polarization [6], thus the lack of a consistent description 
of the polarization in the simulation of the quarkonium 
may result in one of the largest systematic errors that 
may affect the measurement precision. Besides, polar- 



ization observables can provide more information on the 
production and decay mechanisms than the unpolarized 
observables. Given these considerations, the polarization 
issue becomes one of the important topics in quarkonium 
physics for theorists and experimentalists. 

Experimentally, the analysis of direct J/ip production 
at hadron colliders is lacking. The measured prompt 
J/ip data include direct poduction and feed-down con- 
tributions from the Xc and ip' by the radiative decays 
Xc -> J/ipl, ip' -> J/ipJ, and ip 1 -> Xcl -> J/ipll- 
Therefore, in order to compare the theoretical results 
with experiment data, one should also investigate the 
Xc and ip' yields besides the J/ip. Moreover, the Xc 
meson has its own phenomenological interest. The ra- 
tio of differential cross sections of Xci to Xci inclusive 
production at the Tevatron has been measured by the 
CDF Collaboration^]. Their results show that the ra- 
tio disagrees with the spin symmetry perspective from 
the leading-order (LO) computation. After including the 
next-to-leading (NLO) QCD radiative correction [8], the 
asymptotic behavior of gp- changes from p^ 6 at LO to 

p^, 4 at NLO for the 3 p| 8 ' channel and hence becomes com- 
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parable with the contribution of S[ 1 at large px- The 
result provides an opportunity to resolve the contradic- 
tion between the experiment and the theoretical predic- 
tion. Recently, the LHCb Collaboration[9] shows that 
their result of ^ CTxc2 stays within the error bar of the 

NLO NRQCD prediction. Surely, as in the J/ip case, the 
investigation of polarization of Xc will also be very helpful 
in understanding the production mechanism in QCD. 

We now briefly review the recent progress of heavy 
quarkonium in hadroproduction. In Rcf.[10], it is found 
that the NLO corrections for the yields of hard part of 
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J/tp in the S[ channel can be as large as two orders 
of magnitude as compared with LO at large pr, while 
NLO corrections for the CO S-wave are small [11]. For 
the P-wavc, the NLO corrections for 3 P^ channel [8] and 

3 pf^ channel[12-14] are found to be very large but neg- 
ative. As for polarization, the NLO correction for the 
CS channel[15] of the direct J/tp alters from transverse 
polarization in HX at LO to longitudinal polarization 
at NLO, which can be understood in collinear factoriza- 
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tion up to the NLO power in [16]. However, even 

after including NNLO* [17], in which only tree-level di- 
agrams at a| are considered and infrared cutoffs Sij tTn i n 
are imposed to avoid soft and collinear divergences, theo- 
retical predictions of CS contributions to the production 
rates and polarizations are still in disagreement with the 
CDF data[4, 5]. Recently, two groups[18, 19] have pre- 
sented their NLO results for the J/tp direct production 
at hadron colliders, but drawn very different conclusions 
due to different treatments for the fit procedure of avail- 
able data, and extracted different CO long-distance ma- 
trix elements (LDMEs). The former[18] uses a global 
fit, while the latter [19] just concentrates on hadropro- 
duction. The LDMEs in Ref.[19] are more satisfactory 
for the hadron collider data both in yields and polariza- 
tion of the J I tp Especially, the predictions of yields at the 
LHC are in good agreement with data[20, 21] up to 70 
GeV, and the polarizations of jpsi at the Tevatron and 
the LHC are predicted to be mainly unpolarized[19]. 

The study of angular distributions of J / tp from the \c 
decay has been done in Ref . [23] and the authors have also 
calculated the spin density matrix elements (SDMEs) of 
Xd and Xc2 at the Tevatron Run I. The angular distri- 
bution formulas have also been derived in Ref. [24]. In 
the present paper, we re-derive the same expressions for 
the J/tp angular distribution in a different way, and we 
further consider the subsequent decay of the J/tp into a 
muon pair. Measurement of the muon angular distribu- 
tion in the rest frame of J/tp is an additional informa- 
tive result for this process. For the presentation of our 
results with phcnomenological application, we have com- 
puted the inclusive Xci and Xc2 production at the LHC 
with the center-of-mass energy 8 TeV. 

The organization of the paper is as follows. In section 
II, we introduce the basic kinematics and conventions 
used in the paper. In the next three sections, we derive 
the angular distributions of the J/tp and the ji + from 
the decays of Xci and Xc2- In section VI, a re- weighting 
method is proposed to extract all of the SDMEs of the 
Xc2- In section VII, we generalize the rotation- invariant 
relations from the vector boson to arbitrary integer-spin 
particles, and we show that these observables are frame- 
independent. In section VIII, we demonstrate the uncer- 
tainty of J/tp polarization from the feed-down of Xc and 
tp' decays. In section IX, we present the numerical results 
and discuss the corresponding phenomenology. Finally, 
the conclusion is drawn in the last section. The general 
expressions of the angular distributions for the spin-1 and 



spin-2 bosons are listed in the appendices. 

II. KINEMATICS AND CONVENTIONS 

In this section, we introduce the conventions and kine- 
matics for our derivations performed in the following sec- 
tions. We consider the decay process of Xcj — > J/fpJ and 
the subsequent one J/tp — > /i + fi~ . Experimentally, one 
is able to reconstruct J ftp's momentum from its decay 
products /J, + fi~ , which can be identified efficiently at de- 
tectors. One can always choose an arbitrary axis as the 
spin-quantization axis ~£ in its rest frame. Generally, the 
polarization vectors for a massive spin-1 particle are 

£q = (\Ji\,Esm9cos(p,Esiii9sm(p,Ecos9)/m, 

e =P7 

e± = — ^(0,=Fcos#cos(/> + isin</>, 
V 2 

=Fcos#sin0 — icos#,±sin6l), (1) 

where 9, (p are the polar and azimuthal angles of s and 
the symbol 7 can be chosen as an arbitrary real number. 
We set 7 = —(p here. For a spin-2 tensor particle, its 
spin wave functions can be constructed from the spin-1 
polarization 4-vectors as 

ef = £ (l,A i; l,A 2 |2,A)e^e^ 2 , (2) 

Ai,A 2 =-l 

where (1, Ai; 1, A 2 |2, A) are the Clebsch-Gordan coeffi- 
cients. Thus we have identities p M e^ = V\i^\ — Vv e> \ — 
(e A )£ = 0andef = er- 

In Xd — > 4 M + M 7: °nc can work out the angu- 

lar distributions of J/ip in a form of the XcJ production 
SDMEs pw and decay SDMEs D xx 

J 

W(9,4>) = ]T pav-Daa'(M), (3) 

A,A' = -J 

where 9, (p are the polar and azimuthal angles of the J/tp 
or the /i + in the Xcj re st frame. Here p\y and D\y 
represent the production and decay amplitudes of the 
Xcj with spin A multiplied by the corresponding complex 
conjugate amplitudes with spin A'. 

Several polarization frame definitions have been used 
in the literature to fully describe the polarization of 
heavy quarkonium. They are the so-called HX (recoil 
or s-channel helicity frame), the Collins-Sopcr frame, the 
Gottfried- Jackson frame and the target frame [25]. In 
HX, is chosen as the flight direction of the quarkonium, 
while in the Collins-Soper frame s bisects the angle be- 
tween pi and — p~2 in the rest frame of the quarkonium, 
where pi and P2 denote the momenta of the two initial 
state colliding particles in the quarkonium rest frame. In 

the Gottfried-Jackson frame ~£ = &r and in the target 
_^ ml b 

frame s = — t=|t . Here the quarkonium means Xc in the 
case of hadroproduction process pp(j>) — > Xc + X. 
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III. ANGULAR DISTRIBUTIONS IN X a -> J/V>7 with 



The general vertex function for a vector or axial vector 
decays into two vectors can be expressed as 

M(V -> V{V 2 ) = ae Va - e* Vl e V2 ■ (-p Vo - p Vl ) (4) 

+ b e v Q ■ e*v 2 e* Vl ■ (Pv a + Pv 2 ) 

+ ce Vi' e v 2 e v ■ (PVi - Pv 2 ) 

+ d e Va ■ (p Vl - pv 2 )e* Vl ■ (pv„ + Pv 2 ) 

e V 2 ' (PVo + PV 1 ) + f i£e Vo e* Vi e^pv 
+ 9 4i • (PV +Pv 2 )i£e Va c' V2 Pv Q Pv 2 

+ h ey 2 ■ (p Vo +p Vl )ie e 



where 



3p ,o - N Xcl 



2Rpi,- 



37V V 



3 N xd - Po,o ' "-"xci 

3 ^Xd - Po,o 
29fpi,_ 



Po,o 



37V- 



< A/1 



\/2(3pi,0 + %-l,o) 



Xcl 



Po,o 



37V- 



Xcl 



Po.o 



,(7) 



where p it j arc SDMEs for the Xcl yields and A^ Xcl 
Pi,i + Po,o + P-1,-1- 



Pv — Pvi + Pv 2 , ep^pfe = £ nvpcrk a 



Specifically, in the case of Xci(l ++ ) - > J/ip{l )> 
the a, b, c, d terms and the h term can be dropped, be- 
cause of parity conservation in QED and no longitudinal 
polarization for the photon. If we just consider the elec- 
tric dipole (El) transition, which is the dominant contri- 
bution according to the velocity scaling rule in NRQCD, 
the g term can also be negligible. Therefore, we calcu- 



late the helicity amplitudes M, 
(with "/" suppressed) as 



for Xd -> J 1^1 



IV. ANGULAR DISTRIBUTIONS IN X c2 -> J/V>7 



M+++ 

M+Q + 

M+a- 
M 0++ 

M00+ 



= m: 



m x 1 e ^ sin 9 
m v , e 3i * sin 



72 



--M-o- 



-M-o+ = 



Kd + m2 //^) sm2 f 



2m 



COS 



2 e 



2m 



- MSo- = 



J/V 

= -m Xcl e- 2 ^cos6 
K.i+ m J/*) e "^ 



2\/2r 



Similarly, for the Xc2 we can write down the general 
vertex function for a spin-2 tensor particle T 2 decays into 
two vectors 



M(T 2 -> V 1 V 2 ) = a e* V2 ■ e T2 ■ e* Vi 

+ b (p Vl -Pv 2 ) ■ e T2 • e Vl € v 2 ■ {-Pt 2 -PvJ 

+ c (p Vl -Pv 2 ) ■ £t 2 • e V2 e Vl ■ (pt 2 +Pv 2 ) 

+ d (p Vl - pv 2 ) ■ e T2 ■ (p Vl - Pv 2 )e* Vl ■ e v 2 

+ f {Pv l -Pv 2 )- £t 2 • (pvj - Pv 2 ) 
( 5 ) ' (Pr 2 + Pv 2 Yy 2 ■ (pt 2 + PVr ) 

+ Levi-Civita terms, (8) 



where a common factor m 2 ci — m 2 ^ has also been sup- 
pressed. The decay SDMEs are obtained by D\\> = 
Ea. w ,a 7 >'aa. w a 7 ^v Aj/+Ai - Using these ingredients, 
we can work out the general form of angular distribution 
of J/ip from the \d decay 



W 



(M) 



N. 



(l + Xg COS 2 



3 + A 

+ A^ sin 2 9 cos 2<f) + Xq^ sin 29 cos <fi 
+ sin 2 9 sin 2<j> + \ S(t> sin 29 sin 4> 



where px 2 — Pv% + Pv 2 ■ Due to parity conserva- 
tion, we drop the Levi-Civita terms in Xc2(2 ++ ) — > 
J/ip(l )7(1 ). b,c,d,f terms can also be ignored in 
consideration of the fact that we only include the leading 
(6) order contribution, i.e. the El transition, and these terms 



are suppressed by (w 2 c2 — m j/^,) 2 as compared to the a 
term, and some terms should vanish when the photon 
is transversely polarized. Thus, the helicity amplitudes 
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M\ Xc2 \ J/4> \ y for \c2 -> J 7V7 become 



M 2 ++ 
M 2 +- 
M 2 — 
M 2 -+ 

M 20 + 
M 20 - 

Mi++ 
M 1+ - 

M w + 



sin 



M-2— = 



MU-+=e 2t4, cos id -, 
e 4 ^ sin 2 9 



= e 2 ^sin 4tf , 
2 



-M_ 2Q _ = -j= — e ? snr - cos -, 

v2m Xc2 m J/v , 2 2 

-^ 20+ = - m > +m ^ e^cos3^int 
V2m Xc2 m J/v , 2 2 

e"^ sin 20 
-M-i- = ^ , 

= 2e i *cos 3 ^sin^, 



-AT 1++ = 



e 3 ^ sin 29 



-M\, = -2e l</, sin 3 ^cos^, 
2 2 

9 . 9 



= <i±^ si „^ (1+2cos9) , 

2v2rav.„mj/,i, ^ 



.M10- 

A4 ++ = Mo— = 
M00+ 



. . Xc2 J/j> 

•M-10- = —7= 

2V2m Xc2 m J/v , 

"-^^ e 2 ^cos 2 g(2cos0- 
2\/2m„ _ m 7 /./. 2 



67V. 



with 

A^ 

\-L _ 



xci - 9(pia + P-1,-1) - 12po,o 
6 ^xc2 + 3 (Pi,i + P-1,-1) + 4 Po,o ' 
2^(8^/02,0 + Sftp-2,0) + gggi^-i 
67V Xc2 + 3(pi,i + P-1,-1) + 4^,0 ' 

6(K/o 2 ,i - Sftp-2,-1) + V6(Sftpi,o - Sftp-1,0) 
6JV Xea +3(pi,i+p_i,_i)+4po,o 
2\/6(Sp2,o - Sp-2,0) + 63pi -1 
6 ^xc2 + 3 (Pi,i + P-1,-1) + 4 Po,o ' 
6(Qp2,i + Sp-2,-1) + \/6(5pi )0 + Sp-1,0) 
6^x02 + 3(pi,i + P-1,-1) + 4p ,o 



= P2,2 + Pi, 1 + P0,0 + P-1,-1 + P-2,-2- 



(11) 



^* _ Xc2 J/ip 

10+ 2v / 2m Xc2 m 7/v , 

! ^(1 + 3 cos 20) 



4^ 



3 sin 2 9 

-M So = _^«2+^) e ^ sin2g , 
8to Xc2 toj/^ 



The angular distribution of J/ip from the Xc2 decay is 



If one wants to include not only El, but also mag- 
netic quadrupole (M2) and electric octupole (E3) con- 
tributions, the angular distribution becomes much more 
complicated. Nevertheless, one can calculate these con- 
tributions easily by just keeping the coefficients in the 
Xd and Xc2 vertex functions independent and following 
the above procedures. After normalization, there are one 
and two remaining coefficients for \d and Xc2 respec- 
tively that should be extracted from experimental data. 
We have presented them in the appendices. However, 
one should be reminded that when the coefficients of M2 
and E3 contributions are not very small, changes of angu- 
lar distributions due to the higher-order multipole terms 
may also be significant [24] and then one should include 
these contributions in the calculations. The unknown 
higher-order multipole coefficients need to be extracted 
from experiment. However, present measurements of 

l)i these higher-order multipole amplitudes are still not suffi- 
ciently precise and even inconsistent [26-29], and includ- 
ing the M2 and E3 transitions also makes the physical 
analysis more complicated but still does not change the 
underlying physics. Hence, only El transition is consid- 
ered throughout the main context, and we put the for- 
mula including the higher-order multipole effects in the 
appendices, which can be used in practical calculations. 

(9) " 

Our results Eq.(7) and Eq.(ll) are exactly the same 
as those given in Refs.[23, 24]. Actually, our derivations, 
which are based on the general effective decay vertex 
functions, are equivalent to those by using the angular 
momentum analysis as done in Refs.[23, 24], since the ef- 
fective amplitudes written by us are also originated from 
general considerations of spins of involved particles. 



N, 



(l + A 



g COS 



3 + A, 

sin 2 9 cos 2<p + Xg^ sin 29 cos 
Xj, sin 2 9 sin 20 + Xg± sin 29 sin 



(10) 



V. SPIN CORRELATIONS IN 

XcJ -> J/ipl -> P + P~7 

We are now in a position to investigate the /i + angular 
distributions from the cascade decay XcJ — > J/ipl — > 
p+p~7- 

We start from a general formalism to study it. For 
the polarization of XcJ, we suppose its quantization axis 
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as st, while that for J/ip is denoted as s|. Then, the 
general form for the angular distribution of the /i + is 



J/il>l i 



(12) 



J/ip-y 



are production SDMEs of the XcJ, 

,-Ms^f~^ l \ M arc the decay amplitudes 1 
with the lower indices representing their spin or helic- 
ity and 0,<j> denotes as the J/ip's polar and azimuthal 
angles respect to it in the rest frame of the \cj an d 
6' , <j>' denotes as the /x + 's polar and azimuthal angles re- 
spect to i| in the rest frame of the J/ip. Convention- 
ally, the indices that appear twice imply a summation 
here, i.e. J Zl J' z — ± J, ±( J — 1), . . . , 0, s Zl s' z — ±1,0 
and A 7 , A M +, A^- = ±. We emphasize that J z , J' z are re- 
spected to it, while s z ,s' z to it- Again, one can also 
use the general angular momentum relations with proper 
rotations to the specific reference axis to calculate the 
angular distribution function W XcJ ^ J ^^^ + ^~ J (e' , <j>'). 
Nevertheless, they are equivalent for reasons that have 
been given in the previous section. 

Now, we consider the first scenario in which it is de- 
fined as the flight direction of the J/ ip in the rest frame 
of \cj. With the same parametrization as in Eqs.(6,10), 
there are five observables in terms of the XcJ SDMEs 



ft". 






\p = 


1 \ Xcl — \ ±Xcl — n 
~3 ~ <l>' ~ ' 


(13) 


\Xcl _ 

A 8><t>> ~ 






\J-Xd _ 


^(3(P^) + S(P X -to)) 








\y = 


1 

13' 




= 


7>/6(R^ + Rp^ 2 ) + 123^ 








\Xc2 _ 

A 8'<t>' — 


V6(Vtp x _f fi - ^p x $) + 24(R^_ 1 




156iV Xc2 




K> c2 = 


7x/6(S^f - 3p£! 2 ) - 129^1! 








\J-Xc2 _ 

A 0'<t>> - 


VE(Zp x _f fi + 3 P f; 2 ) + 24(^-_ 1 


+ 3p*f) 


156iV Xc2 


5 


with N XcJ 


= Ea=-j Pxx '■ In thc language 


of angular 



momentum theory, the spin correlations between Pj c /ji 



with quantization axis it and p XcJ s ? J ^ J with quantiza- 
tion axis it can be written as 



Jhh = ^ / <m['. ^^IV (14) 

(l,A 7 ;l,s z |J, Ji z )(J, J 2z |l,A 7 ;l,4)Br(xcj -> J/V>7). 



where the summation is also implicitly understood with 
A 7 = ±, J z ,J' z ,Ju,J2z = ±J, ...,0, and Br(x c j -> 
J/ipj) represents the branching ratio of XcJ decays into 
a J/-0 plus a photon, and T>j jr = T>j j, (—(p,0,<p) — 

ei^Jz-J'^dj^j,^), with d J j^ Jt (e) are the well-known 
Wigner d-function 



m'm(J+j' z , J— J z ) 

*Lj>S )= E (-) k - J '- +J - 

k— maxfO, J' z — J z ) 

^{j + j> z )\(j-J' z )\{j + j z )\(j-j z )\ 

(J + J' z - k)\k\{J -J z - k)\{k -J' z + J z )\ 

2J-2k+j' z -J z / n\ 2k-J's+J z 



(15) 



X cos 



One can easily verify that after substituting p x °^ J ^ 1 
calculated from the above equation into the well-known 
angular distribution of the muon from J/ip — > p + pT 



^•//^ + Po,o 

N J/4> + Po,o ' 
V2(Kp{^-Kp%) 

N J/iP + Pofl 



(16) 



the expressions in Eq.(13) are recovered. 

In the second scenario, it is chosen as the same as it 2 . 
The spin eigenstate |1, s z ) with respect to it for the J/ip 
is no longer its helicity state like in the first scenario. 



1 Note, we use general vector current amplitudes for the J/ip de- 
cays into a muon pair. 



2 Note, it is just the case of the J/ip feed-down from the \cj decay. 
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With a direct calculation following Eq.(12), we obtain 



An, 



\J-Xcl 

V 

\J-Xcl 



37V V 



n Xcl ' 
Xcl 



(17) 



Xcl ■ 



3^xd 



3N- 



Xcl 



Xcl 

Po,o 



6A V 



„Xc2 ' 



Xc2 



\Xc2 _ 



\J-Xc2 _ 
A ,h> — 



\J-Xc2 _ 



Xc2 ' 



67V Xc2 +3(p? i f +^f i _ 1 ) + 4^ 

6A r xc2+3(p^+ P ^_ 1 )+4p^ 



6iV- 



Xc2 ^ 
CV„Xc2 



p x t_i: 



■3p*?,-2) + v^(3/#?+^J 



6A r X c2+3(p?; 1 2 +P^t_ 1 )+4p^ 
One can also derive them by combining 

£;£ J /^« ^ £<M a ;M,|J,J a ) 

J* =±1,0 J^J^ 

(J,^|l,Z 2 ;l,s>^Br( XcJ J/Vl), (18) 



and Eq.(16). The expressions in the scenario are also 
exactly the same as those in Ref.[24]. It is interesting 
to see that these observables are mildly corrected by the 
additional M2 and E3 transitions. 

Similarly, one is able to use the same procedure to 
derive the formulas for other scenarios. 



VI. DETERMINATION OF p JziJ , (\J Z - J' z \ > 2) 
FOR X c2 WITH REWEIGHTING METHOD 

From Eqs.(ll, 13,17), we find that the SDMEs 
Pj z ,J'(\Jz — J' z \ > 2) of Xc2 cannot be measured from 
these angular distributions. The fact that the coefficients 
of these SDMEs are suppressed by v 2 or (m Xc2 — mj/^,) 
makes the measurement of these polarization observ- 
ables difficult. In this section, we propose a reweighting 
method to measure these SDMEs. 

From Eq.(9), one can see that this fact is the origin 
for the cancelation of the transverse J/tp and the longi- 
tudinal J/ip from the \c2 decay If the probabilities of 



the transverse and the longitudinal parts are made dif- 
ferent by reweighting, we will have opportunity to split 
pj z ,j> (| J z — J' z \ > 2) without any suppressions. The only 
tradeoff is that one should also measure the polar angle 
0' of the fi + from the J/ip subsequent decay 3 . As is well 
known, the polar angular distribution of the p + from the 
transverse polarized J/ip decay is ( 1+c ° s e ) ; while that 
from the longitudinal one is (1 — cos 2 9'). Therefore, af- 
ter integrating over solid angles of the p + , the weights 
for the transverse and longitudinal J/ip are equal, i.e. 

Dj„J> = EA J/+ =±,0,A,=±^J-.A J /*,A 7 ^J i ,A J ^,A,. and 

it then results in the cancelation of those two parts in 
Pj zi j'(\J z — J' z \ > 2) of \c2- In order to avoid these 
cancelations, a direct way is to measure the 9' of the 
p + and then assign every reconstructed p, + event from 
Xc2 — > J/ipJ — > P + M~7 an extra weight w(cos 2 9'). For 
example, if one chooses w(cos 2 9') oc cos 2 6', because of 



/ 
/ 



d cos <9' (1 - cos 2 6>')w(cos 2 0') 
, (1 + cos 2 9') 



dcosfl' 



the decay SDMEs become 



w(cos 2 6') = 1 : 2, 



D Jz ,j, = 2 Yl M ^ 

+ ■ M ^.0,A 7 A4}, i o,a 7 



(19) 



A-,=± 



and the cancelations do not occur. Without loosing gen- 
eralization, we assume that 



/ 



dcos<9' (1 - cos 2 6»')w(cos 2 0') 

(1 + cos 2 9') . 2nl . 
dcostf' -w(cos 2 6>') = r L : r T , 



and 



D J Z ,J' Z = r T Mj *' 



KA* 

Aj/4,,a~, jvi J' z ,\j/^,,\-, 



^J/ip— ±>A 7 — ± 



+ r L J2 M ' "A M 'r,.. 



(20) 



A 7 =± 



Hence, the angular distribution of J/ip from the Xc2 de- 



3 Note that we are permitted to assume the z axis is in the flight 
direction of J/ip in the rest frame of \c2- 
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cay is J' z \ > 2) if rx 7^ Tl- It is also interesting to note that 

using the re- weighting method, we are not only able to 

yyXc2^J/^7(g) ; 0) oc 1 + \ e cos 2 9 + \ 29 cos 4 9 extract p Jz ^ z (\J Z - J' z \ > 2), but also have the chance to 

+ A sin 2 9 cos 20 + Xj sin 2 9 sin 20 determine all of the SDMEs, which is not met in Eq.(ll). 

For the measurement of the muon angular distribution 

+ A 00 sin 29 cos + A e0 sin 29 sin from thfi ^ decay; the similar wgight w(cog2 ^ can bc 

+ A 2 sin 4 9 cos 20 + A 4 ^ sin 4 9 sin 20 assigned to every event in the rest frame of the \ c2 . 

in 29 -in 2 9 c s 6 Combining the measurements with different weights, one 

2e ^ may also have opportunity to extract all of the produc- 

+ \ 2e4> sin 29 sin 2 9 sin tion SDMEs of Xc2 in the investigation of the produc- 

+ A300 sin 29 sin 2 9 cos 30 il0n mechanism. At the end of this section, we want to 

2 point out that the rewcighting method is also useful in 

+ A304, sin 20 sin 0sm30 t ^ c cx t ra ction of the coefficients in Eq.(21) on the ex- 

+ A40 sin 4 9 cos 40 perimental side when one does an 1-dimcnsional bin fit- 

, 1 . 4 „ . , , s ting rather than the 2-dimcnsional bin fitting. A weight 

+ A4, sm 4 flsm40. (21) , b 2 a\ 1 1 1 a\ 2 a\ u r a <- 

4 * > r v ' w{cos 9) (e.g. w{cos 0) = cos 0) can be applied to sep- 

Thc explicit expressions for the coefficients are arate A< ^' A ^ and Xw ^- 

N Xc 2 — P2,2 + Pl,l + P0,0 + P-1,-1 + P-2,-2, 

i? = 3(rr + r L )N Xc2 + 3r T (pi,i + P-1,-1) + (7rr - 3rL)p ,o, 
6r T N Xc2 - 9r L (p 1A + P-1,-1) - 6(5r T - 3rz,)p ,o 

A " = ' 

A 2e = (rr _ rL) 3^ 2 -15(pi.i+P-i.-i) + 15po,0 | 



A = [2\/6(4rT-3rz,)(Kp2,o + »P-2,o) 

-6(2r T - 3rr,)Kpi,-i] /fl, 
A^ = - [2^(4^-3^)^2,0 -3p-2,o) 

-6(2r T - 3ri)9pi,-i] /.R, 
\ 04 , = [6(2r T -r-i)(»p 2 ,i -»p- 2 ,-i) 

-\/6(2r T - 3ri)(»pi, - Kp-1,0)] /-R, 
A 4 ^ = - [6(2rr - rii)(3p2,i + 3p-2,-i) 

-\/6(2rT - 3r L )(9pi,o + 3p-i,o)] /-R, 

. . . 24»pi,_i - 6i/6(Kp 2 ,o + Kp-2,0) 
A 2 ^ = (r T -r L ) 



R 



-{tt - r L ) 



249pi,_i - 6\/6(9p2,o - 9p-2,o) 



A 2 9^> = 6(r T - ri) [-(Kp 2 ,i - Kp-2,-1) 
+\/6(»pi,o - Kp-1,0)] /i?, 

A29* = 6(r T - r L ) [(Sp 2 ,i + Sp_ 2) -i) 
-\/6(9pi, + 9p-i,o)] 



= 6(r T — r L 



5ftp2,-i - SRp-2.i 



± at s «P2,-1 + 3p-2,l 



R 



(r T - r L ) 



A44, = -(rT-r L ) 



63fy>2,-2 
6$5p2,-2 



7? 



(22) 



From the above equations, we find that one can measure 
A300, A 4 -^, A 4 0, A 4 ^ to determine the values of Pj z ,j z {\Jz~ 



VII. ROTATION INVARIANT RELATIONS FOR 
ARBITRARY SPIN-INTEGER PARTICLES 



The partonic Drell-Yan process in perturbative QCD 
obeys the well-know Lam- Tung identity [30], which states 
that the coefficients Xg and A^ of the lepton angular dis- 
tribution from the Drell-Yan process satisfy Xg+AX^ = 1. 
Its theoretical relevance is that the relation remains un- 
changed up to 0{a 2 s ) corrections[31] and receives rela- 
tively small corrections even by resummation[32]. The 
distinctive feature of the identity is also independent of 
the chosen orientation of the spin axis. Later it was 
pointed out that the rotation invariance of Lam- Tung 
relation is a general consequence of the rotational co- 
variance of J = 1 angular momentum eigenstates by the 
authors of [33]. They presented an expression formally 
analogous to the Lam- Tung identity for a J = 1 boson 
decays into a fermion pair with the only assumption that 
the spin-quantization axis should be set in the produc- 



tion plane, i.e. F\ = 



3+A e 



The new observable Fi 



is rotation invariant in the production plane. The condi- 
tion that the spin-quantization axis is in the production 
plane is indeed fulfilled in HX, the Collins-Soper frame, 
the Gottfried-Jackson frame and the target frame. From 
Ref. [33] , we know that the rotation- invariant property of 
Fi is guaranteed from a relation of the Wigner functions 
^1 mW + ^-1 mW — $\M\,i- In the section, we want to 
generalize the relation to the arbitrary spin-n (n is an 
integer) particles. We straightforwardly write down the 



8 



linear identities for the Wigner functions: 



{k,m-k,m\2k,2m)df mM {e) 

m=—k 

M M 

= {k,^]k,—\2k,M)5 modi M,2)fi, n = 2k, 
k 

^2(2k + l-m, 0; m, 0|2fc + 1, 0) 



m=0 
j2fe+l 



(<?) + 

= (^ + M;^M + M|2 fc + 1,0) 



3 mod(M,2),l; 



n = 2fc + l, (23) 



where k is an non-negative integer. The amplitudes with 
respect to a chosen polarization axis can be symbolically 
denoted as \n) = Y^m=-n a m\ n , m ), where \n, m) is an J z 
eigenstates with the eigenvalues m = —n, —n + 1, ...,n, 
and a m is the production amplitude, i.e. Pj z ,j> z = aj^a}/ . 
From Eq.(23), we can immediately draw a conclusion 
that the linear combinations of amplitudes 



b2k = ^ (fc, m; fc, m|2fc, 2m)a 2m , when n = 2k, 
and 



m— — fc 



6 2fe+1 = ^ (2fc + 1 - m, 0; m, 0|2fc + 1, 0) 



m=0 



(a*2k+l 



■Om-i-2t)i when n = 2k + l, 



are invariant under the rotation in the production plane. 
Therefore, the observables like F n defined as 



m=—n m=—n 

are rotation- invariant, where i?„ is a normalization fac- 
tor to ensure < F n < 1. Or in a more extend 
sense, the functions of F n are rotation- invariant. The 
F n can be expressed into the rational functions of the 
coefficients in the particles' decay products' angular 
distributions^. g.,Ae,A<£, etc). Specifically, for the spin-1 
particles, the observable is 



F 1 = 



1 



\ai + a_i 



2|ai| 2 + |a | 2 + |a_ 



(25) 



while for the spin-2 particles, its expression is 

F 2 = l — 
2 3 |a 2 | 2 + + |a | 2 + |«-i| 2 + |a_ 2 | 2 



\a-2 + V | a o + a^ 2 \ 2 



(26) 



Two concise examples of F\ are the J/ip decays into two 
muons and the \d decays into a J/ip and a photon. For 
the J/ip, 



F[ 



j/i>->-ij. + n~ _ 1 + \q' + 2A^' 



3 + A 9 , 



(27) 



which has been presented in Ref.[33], while for the Xci, 
one can derive 



xci-*-jyv>7 



1 - A e - 4Aj 
3 + Xe 



(28) 



from Eqs.(7,25).The \c2 provides an example of the spin- 
2 particles. In fact, the complete angular distribution 
of the Xc2's decay product J/ip is Eq.(21) in stead of 
Eq. (10). However, the terms absent in Eq.(10) are sup- 
pressed as mentioned above. Hence, the spin information 
in Eq.(ll) is not sufficient. We have proposed a method 
to extract the full spin information of the Xc2 in the last 
section. With Eq.(22) and Eq.(26),one obtains 

tpXc2^J/4>7 n i + n 2^0 + n zX 2 e + n 4 X ct) + n 5 A 2 + n^X^ 
-fS — 



d\ + d 2 Xg + d 3 X 2 e 



m = 

n 2 = 

n 3 = 

n 4 = 

n 5 = 

n 6 = 



3' 

2(8r T + r L ) 
3(10r T -3r L )' 

§r\ - 3r T r L - r\ 
(10r T -3r L )(r T -r L )' 
4(7r T + 3r L ) 
3(10r T -3r L )' 
(2r T -3r L )(7r T + 3r L ) 
3(10r T -3r L )(r T -r L ) : 
7r T + ir L 
3(r T -r L )' 



a 15 a 5 a 3 
di = -,d 2 = -,d 3 = -. 

In particular, if the weight w{x) — x, then r T — 2, r L = 1 
and 

i7Xa2->J/tPi _ 4 2 + 2Ae + 3A 2 e + AX^ + X 2l p + nX^ 



(29) 



15 + 5A e + 3X 29 



(30) 



These frame-invariant relations can be extended to the 
study of other bosons or mesons. We also suggest the 
experimentalists to measure these obsevables to make a 
cross-check of their extractions of the angular distribu- 
tion coefficients in different frames. 



VIII. UNCERTAINTY OF J/ip POLARIZATION 
FROM FEED-DOWN 



The CDF data for the J / ip prompt production include 
not only direct J/ip production but also the feed-down 
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contributions from Xcj and ip'. However, the recent 
NLO calculations of J/tp polarization both in Ref.[19] 
and Ref.[18] are devoid of the feed-down contributions. 
Though the LO NRQCD prediction of the feed-down to 
the J/tp polarization in [34] was found to have a mi- 
nor impact on the final LO result, one may still doubt 
whether the NLO feed-down effect on Xg, 4 of the J/tp 
can be neglected since the NLO correction to the P-wave 
is large [8]. In this section, we will estimate the possible 
uncertainty of the J/ip polarization Xg> arising from the 
feed-down of the Xc and ip' decays. 

The calculation of the prompt J/tp is a little compli- 
cated. In general, prompt data are composed of four 
parts, i.e. the direct production of the J/tp, the single- 
cascade decay of the \c and the ip' and the double cascade 
decay of ip' — ¥ Xcl — ► J '/"077- The direct production 
of the J/tp is straightforward as has been done e.g. in 
Refs.[18, 19]. Relation between the production SDMEs 
of the Xc and the production SDMEs of the J/tp from 
the Xc decay has been given in Eq.(18), while the rela- 
tion between the production SDMEs of the ip' and those 
of the Xc from the ip' decay is 



J>'-*Xc 

E 

=±1,0 



(1, l z] 1, s z \J, J 2 )(l, l z - 1, 4| J, J' z )pf z [31) 



Combine Eq.(31) and Eq.(18), the double-cascade decay 
component can also be calculated. At last, the single- 
cascade decay of the ip' contribution to the J/tp can be 

treated the same as the S[ to the J/tp, because of the 
same selection rule in the electromagnetic dipole transi- 
tion and the chromomagnetic dipole transition. One can 
easily verify it that the above procedure gives the same 
formula as in Ref. [34]. 

Next, we just divide the contribution of the prompt 
J/tp as the direct part p d , which we already know at 
the NLO level and the unknow feed-down part pJ . We 
assume the polarization observables as Xg, and Xg, re- 
spectively. Moreover, the fraction of the J/tp cross sec- 
tion from the feed-down respect to the prompt cross sec- 

- ap k+*toB with p = 



tion is denoted as r, i.e. 



P 



f 



P f o,o 



J- 
Po.o 



and 



1 - AS 



2pi 



^o,o 5 

2pi,l+l»0,0 

po,o- Hence, 



Pi Pia = 



1 + >l 



3 + Xf 



Ptx = (1 - 



1 + 4 
3 + X% 



P(32) 



Note, to be consistent throughout the context, the polarization 
observables of the J/ip or the angular distributions of the muon 
are all denoted by an extra prime. 

Note that we use the symmetry property P^Lx _\i — 
(— ) P\ \i m hadroproduction which is guaranteed by par- 
ity invariance[35]. 
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FIG. 1: Possible impact of the feed-down on the prompt J/tp 
polar asymmetry polarization Xg, with respect to various di- 
rect J/tp polarization X%. Here r = 0.3 is assumed, and the 
curves correspond to J/tp polarizations from the feed-down 

x f e , = ±i,o. 



and 



r | 1 — r 

3+\i, 



(33) 



Two specific cases are analyzed firstly, i.e. Xg, — and 
Xg, = 1, which are the approximated values of the J/tp 
polarization at the Tevatron in HX in Ref. [8] and Ref. [18] 
in pt > lOGeV regime respectively. We obtained the po- 
larization of the prompt J/tp as —^fz < Xg' < zrh. and 



2-6r 



< Xg 



< 1 when varying Xg, from —1 to 1. If we just 
take r = 0.3, which is the usually used value by people, 
the feed-down contribution may change Xg> from 0.24 to 
—0.39 when Xg, = and from 1 to 0.08 when the po- 
larization of the direct J/tp is fully transverse. Hence, 
it seems the impact of the feed-down is more significant 
in polarization prediction in Ref. [18] than in Ref. [19] es- 
pecially when Xg, = — 1. Thus far, the main conclusion 
in Ref. [19] is expected to be not changed too much even 
after including the feed-down contributions as long as 
we use the same LDMEs. In Fig.(l), we establish the 
curves of Xg> with respect to Xg, with Xg, = ±1,0 with 
r = 0.3. The prompt J/tp polarization Xg, can only be 
within the upper and lower curves. Fig. 2 shows the im- 
pact of the feed-down from the Xc and the ip' decay into 
the prompt J/tp polarization Xg, in HX at the Tevatron 
with \/S = 1.96TeV and \y J/4 ,\ < 0.6. The values of 
LDMEs are demonstrated in Tab. I, which is obtained by 
fitting the NLO result with the Tevatron data in Ref. [19]. 
Only central value of the direct J/tp polarization is plot- 
ted. From this figure, we see the feed-down of Xc and 
ip' do not change our result much and our conclusion in 
Ref. [19] is still valid. 



3r 
' 2+r 

from 
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JpJ/W^m^ (oW^Wjj (O^^ 81 )}/™ 2 , The relations between SDMEs of 3 Sf ] and those of ^j 11 

GeV 3 10" 2 GeV 3 10" 2 GeV 3 10~ 2 GeV 3 are similar as in Eq.(31) and are specifically 

1.16 8.9 0.30 0.56 ~ 

3 <?M . v 



TABLE I: CO LDMEs for J/ip from Ref.[19], which is ob- 
tained by fitting the differential cross section and polarization 
of prompt J /if) simultaneously at the Tevatron [5]. The CS 
LDME is calculated by the B-T potential model in Ref.[36]. 



(1, U 1, s a | J, J z )(l, l z ; 1, 4| J, J' z ) Ps ^U) 



direct 

prompt, |A^ ff |-l 




FIG. 2: The J ftp polar asymmetry polarization Xgi in HX 
at the Tevatron with the LDMEs in Tab. I. The dashed lines 
represent the upper and lower limits of the possible feed-down 
contribution to \g> with r — 0.3. CDF data are taken from 
Refs.[4, 5]. 



IX. PHENOMENOLOGY 

In this section, we are now in a position to present our 
results for the inclusive Xd and Xc2 production at the 
LHC with y/S = 8TeV. We also present individual Fock 
states' contributions to SDMEs for the convenience to 
the readers who want to use different LDMEs by using 
our results with simple re-scaling. 

We use our automatic matrix element generator of 
the heavy quarkonium amplitudes to calculate all of 
the SDMEs under various conditions. The generator is 
built up after we rewrote the core codes of the pub- 
lished HELAC[37, 38], PHEGAS[39, 40], RAMBO[41] 
and VEGAS [42]. The program has been tested by a 
number of benchmark processes and has been used in 
the calculation of J/ip polarization successfully [19]. 

The input parameters in our calculations are: 



1) m c = 1.5GeV,m Xc = 2m c = 3GeV. 

2) VS = 8TeV, |y Xc | < 2.4. 

3) PDF set is CTEQ6L1[43]. 



Pt- 



4) fl r = [If 



5) (OX-'CH 11 )) = 3(aJ ti> 2 "° li?', 



\R' P (0)\ 2 



2.2 x 10" 3 GeV J 



4tt 



0.075GeV 5 



with 



• 2,7+1 



Fig. (3), Fig. (4) and Fig. (5) present the transverse mo- 

3 "W 3 PW and 



mentum distributions of 4^ and for 



3 f8l 

S[ respectively in HX, the Collins-Soper, the Gottfried- 

2 , ""^'^ iLJ 



Jackson and the target frames. For P 2 , there is 
contribution too. The cross section in each channel 
( 



dfftot _ 
dp T 



2^ + ^ for 'W and 3 s[ 8] 



' dp T 



dp 



drJtot _ o dfT22 
dpx dpi 



■ i > i 1 fo r *2^) * s a ^ S0 presented, which will be 

used in the normalization of the non-diagonal SDMEs in 
Figs. (6,7,8). We should stress one point that in Fig. (5) 
only (O x -°{ 3 sf ] )) = 2.2 x 10~ 3 GeV 3 are used. In the 
rest of the figures, we demonstrate the angular distribu- 
tion observables derived in Eqs. (7,11,13), where besides 
NRQCD predictions, the CS components are also demon- 
strated. The results show that CO makes a significant 
change compared to the result with CS only. 
Some features should be emphasized: 

1) From the curves of the total cross sections in 
Figs. (3,4,5), we may conclude that the CS dominates in 
the low transverse momentum region, while the CO may 
dominate whenpy increases because gluon fragmentation 
processes [44] become important. The CS pr distribution 
may receive significant contributions from the higher or- 
der radiative corrections [8]. The full NLO predictions for 
Xcj polarizations will be presented later. 

2) From Figs. (3,4,5,6,7,8), we observe that the results 
of pi.j in the Gottfried- Jackson frame and the target 
frame coincide within error bars, aside from a factor 
(— 1)' _J , which is also pointed out in Ref.[23]. For di- 
agonal elements, i.e. pi^, the curves in the Gottfried- 
Jackson and the target frames are similar to those in the 
helicity frame, while many distinctions are found in the 
Collins-Soper frame. Hence for the polarization observ- 
able Xg, the helicity and the Collins-Soper frames are 
enough. Moreover, from Fig. (5), we see that the longi- 

3 [81 

tudinal cross section of the S[ channel is the largest 
in the Collins-Soper frame when px » ^rn c ^ which is 
consistent with the statement in Ref.[45]. 

3) In the helicity and the Collins-Soper frames, the 
SDMEs pi j with \i — j\ odd are almost zero. There- 
fore, they should be measured in the Gottfried- Jackson 
and the target frames as well as Ag^ and Xe'p shown in 
Figs. (11, 14, 15, 17). For p i} j with \i — j\ even, the mea- 
surements in the helicity and the Collins-Soper frames 
are much better than in the Gottfried-Jackson and the 
target frames. 

4) In Figs. (15, 16, 17), we only consider the first scenario 
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in the previous section, i.e. Eq.(13) 6 . In this scenario, 
the El transition just predicts fixed values of Ag, cl , X^f 
and Xgp 2 . On the experimental sides, one can measure 
these observables to see the contributions of suppressed 
transitions like the M2 and E3 parts. These precision 
measurements will improve the test of the validation of 
NRQCD factorization. 

5) The observables defined in Eqs.(25,26) for the Xd 
and the Xci are shown in Fig. (18). The figures show that 
these observables are indeed frame-independent. Hence, 
it would be interesting to measure these observables at 
the LHC. 



X. SUMMARY 

Finally, we draw our conclusion. The upgrade of the 
integrated luminosity at the LHC will not only allow 
us to measure the polarizations of S'-wave quarkonium 
states like J/ip and ip' , but also the angular distribu- 
tions of decay products from the P-wave states Xcj- 
This has opened up new opportunities to further test 
NRQCD factorization and quarkonium production mech- 
anisms in general. In this paper, we present a general 
framework to investigate the polarizations of XcJ and to 
study the impact of these P-wave states to the polariza- 
tion of J/ip by considering the decay vertex functions. 
Specifically, we derive the polar and azimuthal angle dis- 
tributions of the J/ip produced from the Xd and Xci 
decays and the muon produced from the cascade decay 
Xd — s - J /''Pi M + M~7- The coefficients of angular dis- 
tributions, which preserve the spin information of the 
Xcj, are expressed in terms of rational functions of the 
Xcj production SDMEs in Eqs. (7, 11, 13, 17,22). We also 
derive the rotation-invariant relations Eq.(24) for an ar- 
bitrary spin-n particle with n being an integer number 
and apply Eq.(24) to the Xd and Xc2- Moreover, in order 
to illustrate the application of our formulas, we calculate 
the Xd and Xc2 polarizations in HX, the Collins-Soper, 
the Gottfried- Jackson, and the target frames at the LHC. 

Note added: While this paper was prepared, a new 
preprint [46] for polarizations of the prompt J ftp and 
ip' production at the LHC and Tevatron appeared. The 
authors extracted another set of CO LDMEs for the J/ip 
other than those in Refs.[18, 19] by including the Xcj and 
ip' feed-down contributions. Their LDMEs result in two 
combinations of LDMEs that agree with that extracted 
in Refs.[13, 19] and their prompt polarization result is 
around the upper limit in Fig. 2 of this paper. 



Note, in the second scenario, the expressions are exactly the 
same as those for the J/ip angular distribution with only El 
transition is considered as shown in Eqs. (7, 11, 17). Hence, the 
polarization observables defined in Eq.(17) have been shown in 
Figs.(9,12,10,13,ll,14). 
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Appendix A: Angular distributions in the spin-1 
boson decay 



values have been measured in Rcfs.[26, 27, 29]. The nu- 
merical values measured are shown in Tabll. Without 
loosing generalities, the rotation-invariant observable i*\ 
defined in Eq.(25) can be written as 



1 - 3S + (1 - S)X e + 2A 
(l-3<5)(3 + A ) 



(A3) 



In the appendix, we list the general angular distribu- 
tion of the decay products from a vector boson V without 
any symmetry assumptions as 

W y (6>,0) oc 1 + A e cos 2 6» + A sin 2 6lcos(/) 

+ Xecf, sin 29 cos + A 4 ; sin 2 9 sin 20 

+ A 00 sin 29 sin ± 2rjg cos 9 

± 2rjg t j > sin cos + 2r\^ sin sin 0, (Al) 



with 



X e 

Xe<t> 
A 4> 

X-L 

Ve 



(1-35) 
(1-3S) 
(1 - 36) 



(l + S)N v + (l-35)p^ 

(l + 5)N v + (l-3S)pV > 
V2(M,o - Mp\ ) 



(1 + S)N V + (1 - 35)pl a 
(1-35) 



(l + 8)N v + (l-3S)pl ' 

V2(ZpY fi + Zp\ ) 

[ , (l + S)N v + (l-35)pV ' 

Pi,i ~ P-i.-i 
\l + 5)N v + (l-36)pl Q ' 

72(^ + Sftp^ ) 
(l + 5)N v + (l-35)pV ' 

V2(ZpY,o - gp^p) 
a (l + S)N v + (l-35)pl Q ' 



= a 



(A2) 



where N\ 



v 



Pi,i+Po,o+P-u 



1 , and parameters a, S arc 
specified by the identifications of V and its decay prod- 
ucts. In particular, a is induced by the parity violating 
interactions in the decay. In other words, it is non-zero 
only when the decay is not a parity conservative process. 
For the J/ip decays into di-lepton, a — 0, 5 = 0, while for 
the pure El radiative transition Xd ~^ J/ipliOi — 0,(5 = 
i. If one also wants to include the M2 transition in the 



Xd decay, 8 should be changed to 



l+za 1 a 2 



where 



,j=i n J=\ 



represent the El and M2 amplitudes with the 
normalization (af =1 ) 2 + (a^ =1 ) 2 = 1 respectively. Their 



Experiment 


a 2 J=i (10- 2 ) 


CLEO[26] 
Crystal Ball [27] 
E835 [29] 


-6.26 ±0.63 ±0.24 
-0 2 +0 - 8 
0.2 ±3.2 ±0.4 



TABLE II: The M2 normalized amplitude a& for \c 
J/07 measured by different experiments. 



Appendix B: Angular distributions in the spin-2 
boson decay 



For the spin-2 tensor particles, the relevant angular 
distributions are more complicated than that in the pre- 
vious appendix. Here we neglect the details, and only 
write down the angular distribution for a spin-2 tensor 
particle T 2 decay in a general form: 



W T2 (9,<P) oc 1 + A cos 2 9 + A 2 cos 4 9 

± Xgtj, sin 29 cos ± A200 sin 29 sin 2 9 cos 

± A 4 ^ sin 29 sin ± A 4 ^ sin 29 sin 2 9 sin 

± A^ sin 2 9 cos 20 ± A 2 sin 4 9 cos 20 

± Xj sin 2 9 sin 20 ± A 4 ^ sin 4 9 sin 20 

± A300 sin 29 sin 2 9 cos 30 

± A 4 ^ sin 29 sin 2 9 sin 30 

± A 4 sin 4 9 cos 40 + A 4 ^ sin 4 9 sin 40 

± 2770 cos 9 + 2t\iq cos 3 9 

± 27700 sin 9 cos + 2ry 2 00 sin 3 9 cos 

± 27700 sin sin ± 2ry 2 J 00 sin 3 9 sin 

± 2?70 sin 2 9 cos cos 20 ± 2r\^ sin 2 cos 9 sin 20 

± 277300 sin 3 cos 30 

± 2^ sin 3 sin 30, (Bl) 
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where the coefficients in the above expression are 



and 



X g = 6[(1 - 3<J - <fi)JVr a 

- (l-7S + S 1 )(p^ 1 +p\_ 1 ) 

- (3-S Q -75 1 )p^ ]/R, 

X 2e = 0. + 56 - 56i)[Nt !1 -Hp% 1 +p\_ 1 ) 

A e0 = 4[2(l-<5 -2 ( 5 1 )(3fJ^\-3ffp T2 2 _ 1 ) 

- V6(2S - 5i)(Rp^ - Slp\ )]/R, 
A 2e0 = -2(l + 5«5 -5 ( 5 1 )[(3fip^-3fJp T2 2! _ 1 ) 

- V6(3tp^ -n P \ )]/R, 

A^ = 4[-2(l - S - 2<fi)($*p£ + SJp^,.!) 

+ ^(2<5 - ft)(3plf + 9f^ a i, )]/i2, 
A^ = 2(1 + 5ft - 5Ji)[(9f^ + 9fpE! 2 2 ,_i) 

- Ve^ + ^io)]/^ 

A = 4[\/6(l + ( 5o-3 ( 5 1 )(3ffp^ + 3ffp T ^ ) 

- 6(2ft - ft^p^/fl, 

A 20 - -2(l + 5«5 -55 1 )[V6(3fJp 2 ; 2 + 3ffp\ ) 

- mpJ^/R, 

A^ = -4[v^(l + <5o- 350(^0 -$Sp% ) 

- 6(2S -S 1 )Qp^_ 1 ]/R, 

A^ = 2(l + 55 -56 1 )[V6(Zp% ) -Zp% ) ) 

- ASSp^yR, 

ftp? i - 5Rp T2 2 i 
A 3e0 = 2(1 + 5ft -5ft) 1 P 24 , 

So 72 + 3/9 T2 
A^ e0 = -2(1 + 5ft -5ft) ^'-^ / ^ 2 ' 1 , 



A 40 = 2(1 + 55 -5ft) 



AL = -2(1 + 5ft -5ft 



R ' 



% = 2[(2 ai + a 2 )(p% 2 ~ p\_ 2 ) 

- 2(a 1 -a 2 )(p^-p\_ 1 )}/R, 
me = -2(2 ai ~ a 2 )[(pT% - p\_ 2 ) 

- 2{p T l \-p\_ 1 )]/R, 

^ = -A[2(a 1 -a 2 )(ftp^+ftp% ! _ 1 ) 

- V6 ai (ftp^ Q + ftp T \ Q )}/R, 
mH = 2(2a 1 - a2 )[?>{ftpl\+ftp\_ l ) 

- V6(ftpl% + ft P \ a )}/R, 

^ = 4[2(ai - aa)^! - SJp^.-i) 

- ^/6ai(9^-^ li0 )]/fl, 

^ = -2(2ai-a 2 )[3(9^-3pVi) 

- V6(Sp% -Sp\ Q )]/R, 

c fern v^P2,0 — 

?70 = -2V6(2a 1 - a 2 ) — ^ 

^ - 2V6(2ai-a 2 ) — ^ -, 

%60 = -2(2ai-a 2 ) ! - , 



vie* = 2(2ai-a 2 ) 



3p£-i - ^P-2,1 



(B3) 



with 



^T 2 = P^ + Pi\ + Pofl + P-\,-l + P-2,-2. 

i? = (l + 5ft + 3ft)A+ 2 

+ 3(l-35 -5i)(p?] , i+ 

+ (5 - 7<y - 95i)po%. 



(B4) 



(B2) 



In analogy, the parameters ot\ and a 2 are vanishing when 
parity is not violated. In other words, these coefficients 
exhibited in Eq. (B3) are zero in the case of the \c2 decay. 
Other two parameters, ft and ft, can be determined from 
the specific processes one considered. For the %c2 decays 
into a J/ip and a photon, through pure El transition, 
^° = To anc * = irj; wnne after including the higher- 
order multipole amplitudes in the radiative transitions, 
the coefficients ft and ft can be expressed as the poly- 
nomial of El, M2 and E3 amplitudes af =2 , a 2 =2 , a{ =2 

ft = [l + 2a{= 2 (VHa^= 2 +2ai= 2 ) 

+ 4at 2 (at 2 + V5a J 3 = 2 ) + 3(a 3 J = 2 ) 2 ]/10, 

ft = [9 + 6a 1 7=2 (y5a 2 7 = 2 -4a^= 2 ) 

- 4a 2 r=2 (a 2 7 = 2 + 2Vba(= 2 ) + 7(fl3 7=2 ) 2 ]/30.(B5) 

Again, the normalization of (a' 1 /=2 ) 2 + (a 2 =2 ) 2 + (a3 =2 ) 2 = 
1 has been implicitly understood. The numerical values 
extracted from the experiments can also been found in 
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Refs. [26-29] and we summarize them in Tab. III. At last, 
we denote F 2 in terms of the coefficients in Eq.(Bl): 



Hence the coefficients reexpressed as 



F 2 = 

m = 

n 2 = 

n 3 = 

714 = 

n 5 = 

n 6 = 

di = 



n\ + n 2 X g + n 3 X 2e + n 4 A^ + n b X 2lt> + n 6 \ 4c f, 



di + d 2 X g + d 3 X 2 e 



4 - 45 - 35i 
18(2 - 45 - 35i) ' 
2 + 25 - 7fr - A5l + frfr + 3£ 2 
6(1 + 5(5 - 55i)(2 - 4(5 - 3<Si) 1 
1 



3(2-4(50-3^)' 
2S - S 1 



(l + 5(5 -5(5i)(2-4 ( 5o-3<5 1 )' 
1 



1 + 5(5 - 55i ' 

15 A 5 A 3 



(B6) 



Experiment 




a J 3 ='(W--') 


CLEO(Fit 1)[26] 


-9.3 ±1.6 ±0.3 


O(fixed) 


CLEO(Fit 2) [26] 


-7.9 ±1.9 ±0.3 


1.7 ± 1.4 ±0.3 


Crystal Ball [27] 


00 Q+ll-6 
00.0_ 2g 2 


(fixed) 


E760(Fit 1) [28] 


-14 ±6 


(fixed) 


E760(Fit 2) [28] 


-u± 8 7 


+ -t 


E835(Fit 1) [29] 


-9.3t 4 ?±0.6 


(fixed) 


E835(Fit 2) [29] 


-7.6±l A ± 0.9 


2.0t 4 ;|±0.9 



TABLE III 

=2 „J=2 



The 



di- 



al 1 for Xc2 



M2 and E3 normalized amplitudes 
J/ip"f measured by various experiments. 



Appendix C: Dilepton angular distribution in Xc 
decay with multipole effects 



A^ 1 = 

\Xcl _ 

A e>4>> — 

\J-Xcl _ 

xr = 

\Xc2 



3 _ (4=1)2"' V -V ~ U ' 

^K =1 ) 2 (Sft(p?,o 1 )-Sft(^i 1 , )) 
4(3 -(ai=i)*)N Xcl 

4(3-(o^= 1 )2)JV Xel ' 
3(l-ll(a^= 2 ) 2 +9(a^= 2 ) 2 ) 
39 + ll(t^= 2 ) 2 -9(a£= 2 ) 2 ' 
(ai= 2 ) 2 (7^p*f + ^1 2 ) + i mp*^) 

)(^ -^ 2 ) 



(C2) 



2(39 + ll( a r 2 ) 2 -9( a3 J=2 ) 2 )^ 2 
13, 



= [V6(l-y (atr - (al= 2 ' 2 
+ 6(4 - 9(a 2 7 = 2 ) 2 - 4(a£= 2 ) 2 )(H/£ e 2 2 ,-i - Kp*? ) 
/ [4(39 + ll( a2 7 = 2 ) 2 -9(a 3 / = 2 ) 2 )iV Xc2 ], 

(a^ 2 ) 2 (7x/6(3p* 2 2 - 3 Po x ;_! 2 ) - 129^) 



2(39 + ll(a^= 2 ) 2 
13, 



9(« 3 J=2 ) 2 )^; 



Xc2 



A^ 2 

= [V6(l-y (« 2 7=2 ) 2 -(« 3 /=2 ) 2 )(^to + ^f;o 2 ) 

+ 6(4 - 9(a 2 J = 2 ) 2 - 4(a 3 J = 2 ) 2 )(3^_ 1 + 3$*)] 

/ [4(39 + ll( a2 7 = 2 ) 2 -9( a3 ^ 2 ) 2 )iV Xc2 ]. 



Again, we find that some polarization obscrvables are 
trivial and therefore a lot of spin information about the 
Xc disappears. However, the observables in the case can 
be applied to the extraction of multipole amplitudes of Xc 
if these observables can be measured sufficiently precise, 
such as at the electron-positron colliders. 



In the appendix, we consider the di-lepton angular dis- 
tributions of the Xc which has been discussed in section 
IX with only El transition. Here, we try to include the 
higher-order multipole effects in Xc — > J/tp'y — > M + A*~7- 
In the section IX, we have discussed two scenarios to 
describe the angular distributions. Following, we also 
discuss these two scenarios and generalize the formula to 
also account for the higher-order multipole contributions. 

First, the SDMEs for the J/tp from the % c decay should 
be generalized as 

J+l 



^E E I^M::^t^ v . 
1=1 \^=±i j 

(l, A 7 ; 1, s z \J, J\ Z ){J, J 2 z\±, A 7 ; 1, s' z ) 
(2/ + l)(a i / = J ) 2 Br( XcJ ^ J/^y). 



J 



In the second scenario, i.e. the z axis in J/tp — > 
is chosen to be the same as the quantization axis in the 
rest frame of the Xcj, the general relation between the 
SDMEs of the Xcj Pj cJ j> an d those of the J/tp from the 

Xc.j decay p* CJ s 7^ J/ ^ 7 with multipole expansion is modi- 
fied as 



JAA7 * E E E {a J r J ?{i,h;i,s z \j,j*) 

1 = 1 l z =-l Jz,J'^ 

(J,4|/,/ z ;l,4)p^,Br( XcJ J/^ 7 ). (C3) 



(CI) It is clear that the coefficients of the angular distribution 
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of the p + should be written as with 
XT = - N ^ P "\ (04) 



\Xcl 



QN Xc2 - 9(p*f + p*?,^) - 12p x ; 2 + 4(7 - 14(a 2 J = 2 ) 2 - 5(a£= 2 )Xo 1 



\Xc2 
\Xc2 

v 

\Xc2 

A 6'4>' 

Aj-, xc2 



El i2i = [(15-2(a 2 ^ 1 ) 2 )iV Xcl 



A^ 1 = ^=i, - (5-6(«r i ) 2 Ko 1 ]/(5-6(ar 1 ) 2 ), 



i? 2 = [2(21 + 14(a 2 /=2 ) 2 + 5(a^ =2 ) 2 )Ar ; 



Xc2 



i? 2 ' / (7-14(a 2 ^ 2 ) 2 -5(a 3 ^ 2 ) 2 ). 

2V6(Sftp*Q 2 + gp^p) + 

R 2 

6(K^f - Hp^-i) + y6(3?p?; 2 - 3?p*f.o) 



i?2 

2v / 6(Sp^; 2 2 - 3p£! 2 ) - e^p^j 



i?2 

2. 

~ R2 



6(3p* 2 2 + 3p Xc 1 2 2 ) + v / 6(3 ; Po'i 2 +^Po c2 i) Similarly, the feed-down contributions should be modi- 

Afl"'^' 2 = — ~ 1 '~ 2 — 0,-1 , ficd to account for the complete multipole effects. 



[1] N. Brambilla, S. Eidelman, B. Heltsley, R. Vogt, 
G. Bodwin, et al., "Heavy quarkonium: progress, 
puzzles, and opportunities," Eur.Phys.J. C71 (2011) 
1534, 1010.5827. For a very recent review of heavy 
quarkonium production, see also G. Bodwin, 
arXiv:1208.5506. 

[2] G. T. Bodwin, E. Braaten, and G. Lepage, "Rigorous 
QCD analysis of inclusive annihilation and production 
of heavy quarkonium," Phys.Rev. D51 (1995) 
1125-1171, hep-ph/9407339. 

[3] E. Braaten, B. A. Kniehl, and J. Lee, "Polarization of 
prompt J/tp at the Tevatron," Phys.Rev. D62 (2000) 
094005, hep-ph/9911436. 

[4] CDF Collaboration, T. Affolder et al, "Measurement 
of J/tp and tp(2S) polarization in pp collisions at 
yfs = 1.8 TeV," Phys.Rev. Lett. 85 (2000) 2886-2891, 
hep-ex/0004027. 

[5] CDF Collaboration, A. Abulencia et al, 

"Polarization of J/tp and tp2S mesons produced in pp 
collisions at y/s = 1.96-TeV," Phys.Rev. Lett. 99 (2007) 
132001, 0704.0638. 

[6] P. Faccioli, C. Lourenco, J. Seixas, and H. K. Wohri, 
"Towards the experimental clarification of quarkonium 
polarization," Eur.Phys.J. C69 (2010) 657-673, 
1006.2738. 

[7] CDF Collaboration, A. Abulencia et al, 
"Measurement of 

(r Xc2 B(xc2 -> J/ipl)/<? Xcl B{xci -> J I in) in pp collisions 
at y/s = 1.96- TeV," Phys.Rev. Lett. 98 (2007) 232001, 
hep-ex/0703028. 
[8] Y.-Q. Ma, K. Wang, and K.-T. Chao, "QCD radiative 



corrections to XcJ production at hadron colliders," 
Phys.Rev. D83 (2011) 111503, 1002.3987. 
[9] LHCb Collaboration, R. Aaij et al., "Measurement 
of the cross-section ratio a{xc2) /o (Xci) for prompt Xc 
production at yfs = 7 TeV," 1202.1080. 

[10] J. M. Campbell, F. Maltoni, and F. Tramontano, "QCD 
corrections to J/tp and Upsilon production at hadron 
colliders," Phys.Rev. Lett. 98 (2007) 252002, 
hep-ph/0703113. 

[11] B. Gong, X. Q. Li, and J.-X. Wang, "QCD corrections 
to J/tp production via color octet states at Tevatron 
and LHC," Phys.Lett. B673 (2009) 197-200, 0805.4751. 

[12] M. Butenschoen and B. A. Kniehl, "Reconciling J/tp 
production at HERA, RHIC, Tevatron, and LHC with 
NRQCD factorization at next-to-leading order," 
Phys.Rev. Lett. 106 (2011) 022003, 1009.5662. 

[13] Y.-Q. Ma, K. Wang, and K.-T. Chao, "J/psi (psi') 
production at the Tevatron and LHC at 0(a*t; 4 ) in 
nonrelativistic QCD," Phys.Rev. Lett. 106 (2011) 
042002, 1009.3655. 

[14] Y.-Q. Ma, K. Wang, and K.-T. Chao, "A complete NLO 
calculation of the J/tp and tp' production at hadron 
colliders," Phys.Rev. D84 (2011) 114001, 1012.1030. 

[15] B. Gong and J.-X. Wang, "Next-to-leading-order QCD 
corrections to J/tp polarization at Tevatron and 
Large-Hadron-Collider energies," Phys.Rev. Lett. 100 
(2008) 232001, 0802.3727. 

[16] Z.-B. Kang, J.-W. Qiu, and G. Sterman, "Heavy 
quarkonium production and polarization," 
Phys.Rev.Lett. 108 (2012) 102002, 1109.1520. latex, 11 
pages, 4 figures. 



24 



[17] P. Artoisenet, J. M. Campbell, J. Lansberg, F. Maltoni, [32 
and F. Tramontano, "T Production at Fermilab 
Tevatron and LHC Energies," Phys. Rev. Lett. 101 
(2008) 152001, 0806.3282. 

[18] M. Butenschoen and B. A. Kniehl, "J/psi polarization [33 
at Tevatron and LHC: Nonrelativistic-QCD 
factorization at the crossroads," Phys. Rev. Lett 108 
(2012) 172002, 1201.1872. 

[19] K.-T. Chao, Y.-Q. Ma, H.-S. Shao, K. Wang, and Y.-J. [34 
Zhang, "J/ip polarization at hadron colliders in 
nonrelativistic QCD," Phys. Rev. Lett 108 (2012) 
242004, 1201.2675. [35 

[20] ATLAS Collaboration, G. Aad et al, "Measurement 
of the differential cross-sections of inclusive, prompt and 
non-prompt J/psi production in proton-proton collisions [36 
at sqrt(s) = 7 TeV," Nucl. Phys. B850 (2011) 387-444, 
1104.3038. 

[21] CMS Collaboration, S. Chatrchyan et al, "J/psi and [37 

psi(2S) production in pp collisions at sqrt(s) = 7 TeV," 

JHEP 02 (2012) 011, 1111.1557. 
[22] M. Butenschoen and B. A. Kniehl, "World data of J/psi 

production consolidate NRQCD factorization at NLO," [38 

Phys.Rev. D84 (2011) 051501, 1105.0820. 
[23] B. A. Kniehl, G. Kramer, and C. P. Palisoc, "chi(cl) 

and chi(c2) decay angular distributions at the Fermilab [39 

Tevatron," Phys.Rev. D68 (2003) 114002, 

hep-ph/0307386. 
[24] P. Faccioli, C. Lourenco, J. Seixas, and 

H. K. Wohri, "Determination of \ c and Xf> polarizations [40 

from dilepton angular distributions in radiative decays," 

Phys.Rev. D83 (2011) 096001, [arXiv: 1103.4882]. 
[25] M. Beneke, M. Kramer, and M. Vanttinen, "Inelastic 

photoproduction of polarized J / psi," Phys.Rev. D57 [41 

(1998) 4258-4274, hep-ph/9709376. 
[26] CLEO Collaboration, M. Artuso et al, "Higher-order 

multipole amplitudes in charmonium radiative 

transitions," Phys.Rev. D80 (2009) 112003, 0910.0046. [42 
[27] Crystal Ball Collaboration, M. Oreglia et al, 

"Study of the reaction ip' —¥ 77J/W Phys.Rev. D25 

(1982) 2259. [43 
[28] E760 Collaboration, T. A. Armstrong et al, "Study 

of the angular distribution of the reaction 

PP -> Xc2 -> J/ih -> e+e~7," Phys.Rev. D48 (1993) 

3037. [44 
[29] E835 Collaboration, M. Ambrogiani et al, "Study of 

the angular distributions of the reactions 

pp — > Xci, Xc2 -¥ J/ih — > e + e~7," Phys.Rev. D65 [45 
(2002) 052002. 

[30] C. S. Lam and W. -K. Tung, "A Systematic Approach 
To Inclusive Lepton Pair Production In Hadronic 
Collisions," Phys. Rev. D18 (1978) 2447. [46 

[31] C. S. Lam and W. -K. Tung, "A Parton Model Relation 
Sans Qcd Modifications In Lepton Pair Productions," 
Phys.Rev. D21 (1980) 2712 . 



E. L. Berger, J. -W. Qiu and R. A. Rodriguez-Pedraza, 
"Transverse momentum dependence of the angular 
distribution of the Drell-Yan process," Phys.Rev. D76 
(2007) 074006, [arXiv:0708.0578 [hep-ph]]. 
P. Faccioli, C. Lourenco and J. Seixas, 
"Rotation-invariant relations in vector meson decays 
into fermion pairs," Phys. Rev. Lett. 105 (2010) 061601 
, [arXiv: 1005.2601 [hep-ph]]. 

B. A. Kniehl and J. Lee, "Polarized J/ip from chi(cj) 
and yj' decays at the Tevatron," Phys.Rev. D62 (2000) 
114027, hep-ph/0007292. 

S. Chung and T. Trueman, "Positivity Conditions on 
the Spin Density Matrix: A Simple Parametrization," 
Phys.Rev. Dll (1975) 633. 

E. J. Eichten and C. Quigg, "Quarkonium wave 

functions at the origin," Phys.Rev. D52 (1995) 

1726-1728, hep-ph/9503356. 

A. Kanaki and C. G. Papadopoulos, "HELAC: A 

package to compute electroweak helicity amplitudes," 

Comput. Phys. Commun. 132 (2000) 306-315, 

hep-ph/0002082. 

C. Papadopoulos and M. Worek, "HELAC - A Monte 
Carlo generator for multi-jet processes," 
hep-ph/0606320. 

C. G. Papadopoulos, "PHEGAS: A phase space 
generator for automatic cross- section computation," 
Comput. Phys. Commun. 137 (2001) 247-254, 
hep-ph/0007335. 

A. Cafarella, C. G. Papadopoulos, and M. Worek, 
"Helac-Phegas: a generator for all parton level 
processes," Comput. Phys. Commun. 180 (2009) 
1941-1955, 0710.2427. 

R. Kleiss, W. Stirling, and S. Ellis, "A NEW MONTE 
CARLO TREATMENT OF MULTIPARTICLE 
PHASE SPACE AT HIGH-ENERGIES," 
Comput. Phys. Commun. 40 (1986) 359. 
G. Lepage, "A New Algorithm for Adaptive 
Multidimensional Integration," J. Comput. Phys. 27 
(1978) 192. Revised version. 
J. Pumplin, D. Stump, J. Huston, H. Lai, P. M. 
Nadolsky, et al, "New generation of parton 
distributions with uncertainties from global QCD 
analysis," JHEP 0207 (2002) 012, hep-ph/0201195. 
E. Braaten and T. C. Yuan, "Gluon fragmentation into 
P wave heavy quarkonium," Phys.Rev. D50 (1994) 
3176-3180, hep-ph/9403401. 

E. Braaten, D. Kang, J. Lee, and C. Yu, "Optimal 
spin-quantization axes for quarkonium with large 
transverse momentum," Phys.Rev. D79 (1009) 054013, 
[arXiv:0812.3727]. 

B. Gong, L.-P. Wan, J.-X. Wang, H.-F. Zhang, 
"Polarization for Prompt J/rp(ip') production at the 
Tevatron and LHC," , 1205.6682. 



